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STATION KEEPING OF HIGH POWER COMMUNICATION SATELLITES 
by Robert R. Lovell and Thomas A. O'Malley 
Lewis Research Center 

SUMMARY 

Station-keeping requirements were determined for a class of high-power synchro- 
nous equatorial communication satellites characterized by large, Sun-oriented solar 
arrays. The unkept satellite orbit was determined when the combined perturbing forces 
of the Sun, Moon, solar radiation pressure (including Earth shadow effects), radio- 
frequency radiation pressure (from transmitting antennas), and Earth's oblateness and 
triaxiality were acting on the satellite. The effects of radio- frequency radiation pres- 
sure and the Earth's oblateness and shadow are negligible. The remaining forces cause 
three appreciable and separable perturbations in the satellite orbit. 

The Sun and Moon cause a nearly linear (over a 5-year period) increase in the 
inclination of the orbit at a rate of 0. 85° per year. Orbit inclination causes the 
satellite to undergo an apparent daily north-south oscillation with maximum latitude 
equal to the inclination. For minimum propellant consumption, station correction is best 
effected by applying two thrust pulses (12 hr apart) with the center of the pulses occurr- 
ing at the intersections of the desired orbit plane and the actual orbit plane. The veloc- 
ity increment (AV) per year requirement for impulsive thrusting is approximately 
■46 meters per second. 

The Earth's triaxial mass distribution causes a long period, large amplitude oscil- 
lation in the satellite longitude position. The AV per year requirement for impulsive 
thrusting is less than 1.75 meters per second. 

Solar radiation pressure changes the eccentricity of the orbit. An eccentric orbit 
causes the satellite to have an apparent daily east-west oscillation with maximum am- 
plitude in radians equal to twice the eccentricity. For minimum propellant consumption 
station correction for eccentricity is best effected by applying two thrust pulses 12 hours 
apart. These thrust pulses control the line of apsides such that the projection of the 
. Earth-Sun line into the orbit plane is coincident with the Earth-perigee line. The AV 
requirement is proportional to the area-to-mass ratio of the satellite. For an area-to- 
mass ratio of 0. 15 square meters per kilogram, an average satellite reflectivity of 
0. 3, and an allowable longitude error of 0. 15°, the AV per year requirement for im- 
pulsive thrusting is 8 meters per second. 



INTRODUCTION 


The next generation of communication satellites will use large, lightweight, Sun- \ 
oriented solar arrays as a primary power source. The satellites will transmit narrow | 
radio-frequency beams from synchronous equatorial orbit to high gain receiver antennas | 
on the ground. The system engineer must make a trade-off between the gain and pointing| 
requirements of the ground receiver antenna and the station-keeping accuracy of the | 
satellite. This is particularly true for missions using frequencies around 12 gigahertz. | 
For example, at 12 gigahertz, with effective receiving apertures as small as 4 square 1 
feet (0. 37 m^), ground receiver antennas will have a half power beam width of 2°. Be- ) 
cause of antenna mounting uncertainties, it may be mispointed by 1/2°. If the satellite J 
drifts off station by an additional 1/2°, then the received signal strength will be reduced 2 
by a factor of two. If the satellite prime power is fixed Thither the ground antenna point-1 
ing accuracy must be improved or the satellite station-keeping accuracy must be in- 
creased. 

The problem of satellite station keeping is not new. The yearly AV requirements 
of 46 meters per second for north-south and 2 meters per second for east-west station 
keeping are documented in a number of reports. The 46-meter-per-seconu AV re- j 
quirement for north-south or inclination control is based on impulsive corrections re- 
quired to counter the perturbing accelerations of the Sun and Moon. The 2-meter-per- 
second AV requirement for east-west or longitude control is based on impulsive cor- 
rections required to counter accelerations arising from the Earth triaxial mass distri- 
bution. For relatively dense, rigid satellites, station- keeping requirements can be 
specified by considering only the previously mentioned effects. But for high- power 
communication satellites (HPCS’s), this is not enough. 

HPCS’s have two characteristics that complicate the problem of determining 
station-keeping requirements. The first of these characteristics is the reflective prop- 
erties of large Sun-tracking solar arrays. Solar pressure causes an acceleration pro- 
portional to the f .’ea-to-mass ratio of the satellite, and the resultant accelerations , 
change the eccentricity of the satellite orbit. As a result of the large Sun -tracking solar 
array, the area- to- mass ratio is sufficiently high that the station- keeping requirements 
due to solar pressure must be considered. 

Second, HPCS’s behave as flexible bodies. Large roll-out solar arrays are very 
flexible and can tolerate only mild accelerations without structural failure. Impulsive 
accelerations resulting from station corrections initiate solar array panel oscillations 
that interact unfavorably with the fine pointing attitude control system. As a result of 
the flexible structure characteristics, the station-keeping impulses must be distributed 
over a longer period of time. 
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This report presents the results of an analytical study of the requirements neces- 
sary to keep a class of 24-hour synchronous equatorial communication satellites on 
station. The class of satellites considered can be described as Sun-tracking flat plates 
with area- to- mass ratios varying from 0. 05 to 0. 15 square meters per kilogram. This 
includes the HPCS series of satellites, which have a dense central body, some relatively 
small reflector antennas, and large Sun-oriented solar arrays. 

The reader is presented with working curves and equations that show the relations 
between spacecraft parameters, station accuracy, and thruster and propellant require- 
ments. The approach taken is to present the reader with a review of the unkept (or what 
happens if one does nothing) satellite station when the satellite is subjected to varicvs 
perturbing forces. The unkept station is described in terms of the classic orbital ele- 
ments. The perturbing forces considered are the Sun, Moon, solar pressure including 
the Earth’s shadow, Earth’s oblateness and triaxiality, and radio-frequency radiation 
pressure. 

After establishing the unkept station, fundamental concepts of correcting the per- 
turbed orbital elements are examined. Radial, tangential, and out-of-plane thrust 
schemes are examined to determine their applicability to controlling the orbital ele- 
ments. Satellite accelerations and characteristic velocity requirements and hence 
thruster size and propellant weights are given for various r athods of station keeping. 
Finally, a sample problem is presented to demonstrate the use of the design curves for 
a typical HPCS mission. 


SYMBOLS 

Symbols used for special purposes are defined where they occur in the report and 
are not included in the following list. 

o 

A/m area-to-mass ratio of satellite, m /kg 

a semimajor axis 

a^jg thruster acceleration, north-south station keeping 

a . thruster acceleration, east-west station keeping due to solar pressure for the 
S1 th 

i n method 

e eccentricity 

e* maximum allowable eccentricity 

ep peak eccentricity that would occur if initial orbit were circular (e Q = 0) and no 

station keeping were applied 
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i inclination 

k (A/'m)(l + a) 

AL allowable longitude error 

AL e allowable longitude error due to solar pressure effects 

AL t allowable longitude error due to triaxiality effects 

M number of days to complete a station-keeping correction of a specified type 

N number of days between the beginnings of successive station-keeping correc- 

tions oi a specified type 

P unit vector from center of Earth to orbit perigee 

P., unit vector formed by projecting P into equatorial plane 

x 

p duty cycle, equal to thruster on-time per orbit divided by orbit period 

Q unit vector orthogonal to P such that P and Q lie in orbit plane 

r radial distance of the satellite from the center of the Earth 

Ar variation of r from the nominal synchronous radial distance 

Ar g variation of r due to solar pressure effects i 

Ar^ variation of r due to triaxiality effects 

-fi 2 

S solar constant, 4.5x10 kg/(m)(sec) 

sk station keeping 

U unit vector from center of Earth to San ! 

_ _ 1 
Uj unit vector formed by projecting U into equatorial plane 

V nominal satellite velocity, 3075 m/sec i / 

AV C velocity increment needed to make a station-keeping correction 

AVjjg velocity increment per year for north-south station keeping 

AV_ 4 velocity increment per year for east-west station keeping (due to solar pres- 

sure) when using the i method 

AV t velocity increment per year for east-west station keeping (due to triaxiality) 

X, Y, Z inertial coordinate system 

/3 eccentricity ratio e*/e 

r 

y longitude of satellite measured from minor axis of Earth's equatorial section 

y Q desired value of y 
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y value of y at the perigee of orbit just prior to an east-west correction 

Jr 

Ay variation of y from y Q where y = y Q + Ay 

Ay p variation of y from y p where y = y p -r Ay p 

Ay varition of y caused by solar pressure effects 
s 

Ay t variation of v caused by triaxiality effects 
A longitude of Sun measured from X-axis 

A mean angular velocity of Earth's orbit about the Sun A= 1.99x10 rad/sec 

8 angular velocity of Earth’s rotation about its axis = 7.29x10 ® rad/sec 

I u gravitational constant of Earth 

a average reflectivity of satellite 

<p latitude of satellite 

if/ angle through which apsidal line is rotated 

u) longitude of perigee measured from X-axis 


SATELLITE AND ITS ORBIT 

Figure 1 shows three typical satellites of the type considered in this report. These 
satellites have a dense central body, some relatively small reflector antennas, and 
large Sun- oriented solar arrays. They can be characterized as Sun- tracking flat plates 
with area-to-mass ratios varying from 0.05 to 0. 15 square meter per kilogram. 

To facilitate continuous Earth coverage to nonsteerable ground antennas, these sat- 
ellites will operate from synchronous altitude. Satellite position control to within ±0. 2 
in latitude and longitude may be required. Typical missions will last for 5 years. 


NATURAL PERTURBATIONS AND UNKEPT STATION 

Because of the effect of various perturbing forces, a satellite in synchronous equa- 
torial orbit will not remain stationary with respect to the rotating Earth, Perturbing 
forces considered in this report arise from the oblateness and triaxiality of the Earth, 
the Sun and Moon gravitational attraction, and solar radiation pressure. Perturbation 
forces arising from radio-frequency radiation pressure from communication antennas 
are neglected because they were found in appendix A to have negligible effect on the sat- 
ellite orbit. 
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Before the methods and requirements of station keeping are analyzed, the resultant 
behavior of the unkept (no station keeping) satellite station will be determined. 

The position of a satellite with respect to the rotating Earth can be completely de- 
termined by three parameters: r, distance from center of Earth to satellite; y, satel- 
lite longitude measured from the minor axis of the Earth's equatorial section; and ip, 
satellite latitude. Letting r Q be the nominal value of r fora24-hoi.: circular orbit 
and y Q be the initial, value of y, the variations Ar and Ay can be defined by the 
equations 


r = r + Ar 
0 

(l) 

y = r 0 + Ay 

(2) 


The deviations of a satellite from a nominally synchronous equatorial or* .t are given by 
Ar, Ay, and <p. An out-of-plane perturbation creates a nonzero orbit inclination which 
in turn causes nonzero values of <p. An in-plane perturbation causes changes in orbital 
period, eccentricity and orientatioi. of the apsidal line which in turn causes changes in 
Ar and Ay. 


Out-of-Plane Perturbations (Changes In (p) 

From reference 1, the effect of the Sun and Moon is to cause the orbit inclination to 
build sinusoidally to a peak of 14. 7° after 26^ years and then decrease to zero after 
53 years. Over the first 5 years, the inclination increases at an approximately linear 
rate of 0. 85 degree per year. The satellite latitude undergoes a sinusoidal oscillation 
once per orbit with amplitude in radians equal to the instantaneous orbit inclination. 
Figure 2 shows a plot of latitude as a function of time. The effect of the Sun and Moon 
is the only significant out-of-plane perturbation. 


In-Plane Perturbations (Changes In Ar anti Ay) 

The Earth's oblateness and triaxiality, the Sun and Moon, and solar radiation pres- 
sure cause in-plane perturbations. The effect of oblateness is to modify the value of r Q 
from that obtained from a spherical Earth model. The effect of the Sun and Moon is to 
cause small oscillations in both Ar and Ay. (From ref. 2, a conservative upper bound is 
3000 meters for Ar and 0.06° for Ay. ) The effect of tTiaxiality is to cause a longttu- 
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dinal oscillation about the Earth's mionr axis. The effect of solar pressure is to change 
tg 

solar pressure are importantaid are analyzed below 

Triaxiality . - From reference 3 the Earth's equatorial cross section is approxi- 
mately an ellipse with an ellipticity of j, where has the value of -1. 816x10"®. 

The minor axis of the equatorial ellipse passes through 74. 6° east longitude and 105. 4° 
west longitude (see ref. 3). These two longitudes are stable points. Neglecting other per- 
turbations, a satellite placed at either of these longitudes (y Q = 0)will tend to stay there. 
Hence, if a satellite is positioned at any other longitude (y Q * 0), it will undergo a longi- 
tudinal oscillation about the nearest minor axis (see ref. 4). The period of this oscillation 
is greater than 2. 2 years and is a function of y Q . The variation of r also undergoes an. 
oscillation of the same period. The amplitude of this oscillation is a function of y Q . 

When y Q is 90°, the amplitude of the Ar oscillation is a maximum (35 000 m). The 
variation of r and y due to triaxiality will be denoted by Ar t and Ay^. Figure 3 
presents \r t as a function of time with y Q = 45°. Figure 4 presents Ay ( as a function 
of time with y Q = 45°. 

Solar pressure . - The effect of solar pressure is to change the orbit eccentricity 
and orientation of the apsidal line. It has a negligible effect on the - rbit period. The 
induced eccentricity causes a daily longitudinal oscillation with amplitude 2e radians. 

The period of this longitudinal oscillation (24 hr) is much smaller than the period of 
longitudinal oscillation due to triaxiality. Solar pressure also causes a daily oscillation 
in orbit radius with amplitude er Q . Jnlike the other perturbations, the effect of solar 
pressure is dependent on satellite parameters. With the assumptions that the satellite 
is a flat plate whose surface is perpendicular to the Earth-Sun line and that the front side 
thermal radiation is the same as the back side thermal radiation, the perturbing accel- 
eration of the satellite due to solar pressure is 

a = -SkU (3) 


where U is a unit vector from the center of the Earth to the Sun, S is the solar con- 
stant at 1 AU, and 


k = (1 + a) — (4) 

m 

or approximately 1. 3 times the area- to- mass ratio for silicon cell solar arrays (as- 
summing average satellite reflectivity a = 0. 3). 

There is no solar pressure when the satellite is in the Earth's shadow. For 275 days 
of the year, the satellite does not enter shadow at all. Over a 1-year period, the satel- 
lite is in shadow 1 percent of the time. The satellite is in shadow no more than 5 per- 
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cent of the orbit even for the worst case (at the vernal and autumnal equinoxes). The 
inclusion of Earth shadow in computer solutions indicated that the effect of shadow is 
negligible. The effect of shadow is ignored in the following analysis. 

Before continuing with the discussion of solar radiation pressure, refer to figure 3 
for definitions of the quantities U, Uj, P, Pj, X, and co. The XYZ reference cooidi- 
nate system is an inertial system with the origin at the center of the Earth. The X-axis 
is toward the autumnal equinox, the X-Y plane is the equatorial plane, and the Z-axis is 
along the Earth's spin axis. Notice that w is not th r argument of perigee, nut rather 
the longitude of perigee measured from the X-axis. When the initial orbit is circular, e 
as a function of time is given to a good degree of approximation by 

e(t) = e !sin — Xtl (5) 

P| 2 I 


where 


e p is given by 


e =£££= 0.022 k (6) 

p VX 

# 

where k is in square meters per kilogram, V is the nominal satellite velocity, and X 
is the mean angular velocity of the Earth's orbit about the Sun (2v > d/yr), Longitude of 
perigee o>(t) in radians is given approximately by 

w(t) = x + - + A Xt; 0 < t < 1 year (7) 

0 2 2 

Equation (7) implies that the line of apsides rotates uniformly in inertial space, making 
a 180° rotation in a 1-year period. Equations (5) to (7) a r e derived analytically in ap- 
pendix B with the assumption that the ecliptic plane, equatorial plane, and orbit plane 
are the same plane. Solutions for e(t) and cu(t), without the coplanar assumption, were 
obtained on a digital computer by using numerical integration. The computer solutions 
agreed well with analytic solutions. Figure 6 presents the computer solutions for 
e(t)/k for two cases. In the first case (starting at the vornal equinox), the maximum 
eccentricity is slightly less than for the second case (starting at the winter solstice) 
When starting at any other time of the year, the maximum eccentricity will be some- 
where between the two maximums shown in figure 6. Figure 7 presents u>(t) and X(t) 
when starting at the vernal equinox. The apparent discontinuity in w\t) at time equal 
to 1 year is resolved by realizing that the line of apsides is undefined when the eccen- 
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tricity is zero. The relation of u to X plays an important role in choosing station- 
keeping techniques. 

Appendix C gives a detailed discussion of the resulting orbit when the initial orbit is 

1 

noncircular. However, one special case deserves attention here. If e = -* e and 

O ci y 

a = X Q> then e(t) and w(t) are given to a good degree of approximation by 


e(t) = - e 
2 p 

(8) 

-M 

ll 

+ 

o 

<< 

II 

-4— > 

3 

(9) 


Thus, the eccentricity remains constant, and the line of apsides rotates uniformly in 
synchronization with the Earth-Sun line. Equation (9) implies that the orbit perigee is 
Sun-oriented. The analytic derivation of equations (8) and (9), using the coplanar as- 
sumption, is included in appendix B. A 24-hour orbit having an eccentricity and longi- 
tude of perigee as given by equations (8) and (9) will be called a Sun-oriented orbit. The 
importance of the Sun-oriented orbit is that its eccentricity remains constant at a value 
of only half the maximum eccentricity obtained from an initially circular orbit. The 
Sun-oriented orbit is analogous to a Sun- synchronous orbit. In a Sun- synchronous orbit, 
a proper combination of orbital radius and inclination will cause the line of nodes to 
rotate with an angular velocity of 360° per year. In the Sun- oriented orbit, the proper 
combination of eccentricity and initial longitude of perigee causes the line of apsides to 
rotate with an angular velocity of 360° per year. Figure 8 presents plots of e(t)/k ob- 
tained from computer solutions using numerical integration. Figure 9 presents u>(t) and 
, X(t). The computer solutions did not use the coplanar assumption. Both figures show 
f close agreement to equations (8) and (9). 


Summary of Perturbations 

Three significant and distinct orbital motions result from the perturbation forces 
acting on a synchronous equatorial communication satellite. They are summarized as 
follows: 

(1) The Sun and Moon cause the orbit to develop an inclination at the rate of 0. 85° 
per year. The induced inclination causes the satellite latitude to undergo a daily sinus- 
oidal oscillation with amplitude equal to the instantaneous inclination. 

(2) The Earth’s equatorial section is approximately an ellipse with minor axes at 75° 
east and 105° west longitude. A satellite positioned at any other longitude will tend to 


9 



drift toward and oscillate about the nearest minor axis. The period of this oscillation is 
greater than 2.2 years and is a function of the initial longitude. 

(3) Solar pressure changes the orbit eccentricity and longitude of perigee. The in- 
duced eccentricity causes the satellite longitude to undergo a daily longitudinal oscilla- 
tion with amplitude 2e radians. The perturbing acceleration of the satellite due to 
solar pressure is proportional to the area- to- mass ratio. 


STATION KEEPING METHODS AND REQUIREMENTS 

In the previous sections, the orbit perturbations were presented in terms of A r, 

Ay and <p. Changes in Ar, Ay or cp will cause changes in one or more of the orbital 
elements a, e, i, and co. Table I presents a summary of the perturbations, their ef- 
fects on the orbital elements, their effects on Ar, Ay, <p, and some comments. 


TABLE I. - SUMMARY OF PERTURBATIONS 


Perturbation 

Effect on the 
orbital elements 
a, e, i, w 

Effect on Ar, 
Ay, (p 

Comments 

Sun and Moon 

i = (0. 86°/yr)t 
Small oscillations 
in a 

(p - i sin 0 E t 
Small oscillations 
in Ar, Ay 

Station seeping needed to 
control i 

No station keeping needed 
to control a 

Oblateness 



Oblateness modifies r 

0 

Nc station keeping needed 
ay long as orbit is nomi- 
nally equatorial 



Triaxiality 

Long- period 
oscillation 
in a 

Long- period 
oscillation in 
Ar t and Ay^ 

Station keeping needed to 
control a 

Solar radiation 
pressure 

1 

If e„ - 0, 

0 

then 

e = e^ sin I Xtl 

p l 2 | 

10 X + 1 + A Xt 

0 2 2 

If e = i e„ and 
0 2 p 

'•i = A , then 

o O' 

e = i e 
2 p 

u) *= ^t + X Q 

Short- period 
oscillation in 

Ar s- A *s 

Station keeping may be 
needed to control e 
and/or to 
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In discussing station-keeping methods and their associated requirements, it is con- 
venient to use equations that give the time rate of change of orbital parameters as a 
function of the station- keeping acceleration. Integrating these equations over the time 
interval of thrusting will yield the change in the orbital parameters. The station-keeping 
acceleration vector can be lesolved into a component R along the radius vector (meas- 
ured positive away from the Earth), a transverse component T in the instantaneous 
orbital plane (measured positive when in the same direction as the orbital velocity vec- 
tor), and a component W normal to the instantaneous orbital plane (measured positive 
to the north). From reference 5, the equations for the time rate of change of the orbital 
elements are 


da 2e sin 0 


dt 


n/ 




1 - e 


nr 


de _ 1 - e^ sin 6 R | 


dt 


na 


na^e 


a^(l - e^) 


- r 


( 10 ) 


(ID 


-- = - V - 1 ~ e2 P QS — R + t - 2 (l + 1 Vsin 0)T (12) 

dt nae nae \ 1 + e cos 0/ 


di - r cos £ yff (13) 

' dt 2JZ 2 

^ na j 1 - e 

X 

5 

11- where n is the orbital angular velocity, r is the orbital radius, 9 is the true anomaly, 

r and cp is the angle from the ascending node to the instantaneous position of the satellite. 

Retaining only first-order perturbations in e and assuming e cos 8 « 2, equations 
(10) to (13) can be reduced to 

£ 



f 


1 , 



Jf 


da _ 2ae sin 9 R ( 2a T 
dt V V 

de _ sin 8 R ( 2 cos 8 ^ 
dt V V 

d(» _ -cos 8 R | 2 sin 0 T 
dt eV eV 


(14) 

(15) 


(16) 
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di 

dt 


cos <? w 
V 


(17) 


In all of the station-keeping methods to be described, there are two thrusting 
periods per orbit of 12 p hours each where p is the thruster on-time per orbit divided 
by the orbit period (24 hr). The first thrusting period is centered about a position in the 
orbit to be called Aj. The second thrusting period is centered about Ag. In all cases, 
Ag is separated from Aj by an angular distance of 180°. Each station- keeping correc- 
tion is assumed to take place over M consecutive orbits. Thus one station-keeping 
correction consists of 2M thrusting periods and a total thrusting time of Mp days. 

In the following sections, the station-keeping requirements are given in terms of 
the AV per year and the thruster acceleration. 


North -South Station -Keeping Methods and Requirements 
Due to the Sun and Moon 


In the previous sections, it was shown that the Sun and the Moon cause the orbit to 
develop an inclination. To control the effects of the Sun and the Moon, one must control 
the inclination of the orbit. To change the orbit inclination an amount Ai, it is neces- 
sary to thrust normal to the orbital, plane. Figure 10 presents a sketch of this maneuver 
for M = 1 and for a duty cycle close to zero (impulsive thrusting). The correction is 
made by thrusting in the south direction in the vicinity of Aj and thrusting in the north 
direction in the vicinity of An where Aj and Ag are the ascending and defending nodes, 
respectively. The change in inclination Ai as a function of duty cycle p can be found by 
integrating equation ( 17 ). 


Ai = 2M 




W dt 


V 


The result is 


Ai = 


2 AV„ sin ££ 
c 2 


Vpn 


( 18 ) 


( 19 ) 
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gives 


Solving for AV 


AV = V Ai — 25 (20) 

2sinH 

2 

Let aj^g be the north-south acceleration caused by the station-keeping thruster. 
Then a^g = W and the change in inclination for a complete correction (firing near 
each node for M consecutive orbits) is given by 


Ai = 


4Ma NS . pn 
— sin 

«E V 2 


( 21 ) 


Recalling that the Sun and Moon cause i to increase at the rate 0.85 degree per year, Ai 
can be expressed in radians as 


Ai = 


0.85 N 
57.3 365 


( 22 ) 


where N equals the number of days between the beginnings of successive inclination 
corrections. Equations (21) and (22) can be combined to yield 


(2.3xl0 -6 )-^ 

%ec 2 ) sin 2! 


EE 

2 


Using equation (20), the AV per year (AV ns ) is given by 

/J*A x 46 — H— 

\sec) 


AV, 


NS 


2 sin-2E 
2 


(23) 


(24) 


The AV Ng is plotted as a function of p in figure 11. In figure 12, Ma^g is plotted as 
a function of p for N = 1, 7, and 60. 
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East-West Station -Keeping Methods and Requirements Due toTriaxiality 

Triaxiality causes a long period east-west or longitude oscillation. Accompanying 
the longitudinal drift is a change Aa in the semimajor axis. The method used to control 
against triaxiality is based on changing the semimajor axis. 

The semimajor axis of the orbit can be increased by some small amount (Aa > 0) by 
thrusting eastward in the vicinity of both A ^ and Ag. Position Aj can be either the 
apogee or the perigee of the orbit. If Aa is negative, then the direction of the thrust is 
westward. Figure 13 presents a sketch of this maneuver for impulsive thrusting when 
Aa is negative. For impulsive thrusting, equation (14) can be integrated to yield 

Aa = — AV (25) 

V c 

Solving for AV c gives 

AV = — Aa (26) 

c 2a 

For nonimpulsive thrusting, 



Figure 14 presents a plot of Aa against Ay^. The dashed portion of the curve 
represents the motion when station keeping is not used. When station keeping is used, 
the satellite will be initially positioned so that Ay t = - AL t , where AL^ is the maximum 
allowable excursion due to triaxiality. The initial semimajor axis will be a Q n . fsee 
point A, fig. 14). The satellite will drift eastward until Ay t = +AL t (point B, fig. 14), 
and then begin drifting westward. When the satel’ite has drifted to a point where 
Ay t = -AL t (point C, fig. 14), the semimajor axis will be a Q - Aa c . Station keeping is 
then used to increase the semimajor axis an amount 2Aa c , bringing the satellite back to 
point A again. The process is then repeated. The station-keeping acceleration level is 
small enough so that the corrections can be done impulsively. The AV per year (AV^) 
for correcting triaxiality is then only a function of the off-longitude y Q (see ref. 4). 
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AV /— \ = 1 - 75 l sin 2 r 0 l 

\sec/ 


(28) 


Figure 15 presents a plot of AV^. as a function of y . 


East-West Station-Keeping Methods and Requirements Due to Solar Pressure 


* 


In the previous sections, it was shown that solar pressure changes the eccentricity 
and rotates the line of apsides of the orbit. It will be shown in this section that the ef- 
fects of solar pressure can be controlled by continuously thrusting directly against the 
Sun, or by controlling the eccentricity, or by controlling the line of apsides. Four 
methods are presented. Appendix D presents station-keeping techniques for changing 
orbital eccentricity and rotating the line of apsides. Derivations of the formulas for AV 
and acceleration requirements for each method are presented in appendix E. 

Method 1 . - In this method the effect of solar pressure is canceled by continuously 
thrusting toward the Sun. The acceleration of the satellite due to the thrust is equal but 
opposite to the acceleration caused by solar pressure. The AV per year is 


AV sl = 


_ 3Sk;r /4 


2X \3 


(29) 


The acceleration level is 


% 


a sl - sk 


(30) 


JTj 




'f 

t 

t 

f 

.V 


Method 2 . - In this method each time the eccentricity becomes equal to a predeter- 
mined maximum allowable eccentricity e* the orbit is circularized. Appendix C shows 
that an eccentric orbit can be circularized by thrusting either collinearly with the orbit 
velocity vector or collinearly with the orbit i dius vector. Appendix C also shows that 
tangential thrusting requires only one- half as much AV as radial thrusting. For that 
reason, the requirements for this method are given for tangential thrusting only. Let 
AL_ denote the maximum allowable longitude excursion due to solar pressure consist- 
ent with the station accuracy requirement. Since eccentricity causes a daily longitude 
oscillation of amplitude 2e radians, e* = ^ AL g (where AL g is in radians). The param- 
eter /3 is defined as the ratio of the maximum allowable eccentricity to the maximum 
eccentricity that would result from an initially circular orbit (assuming no station keep- 
ing), that is, 


/ 
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(31) 


0 


e* 


Clearly, if 0 > 1, no station keeping is required. If AL g is in degrees and k is in 
m^/kg, 0 can be expressed as a simple function of AL g and k. 


0 = 



(32) 


Figure 16 presents a plot of e/e p , o>, and X as functions of time when method 2 is 
used. The AV per year is 


AV 


_ 3Sk?r 
s2 = 2X 



The acceleration level is 


a s4 = 



(33) 


(34) 


Figure 17 presents a plot of AV^/k as a function of p with 0 as a family parameter. 
Figure 18 presents the corresponding plot for (M/k)a g 2 . 

Method 3 . - In method 3, the line of apsides is rotated each time the eccentricity 
becomes equal to the maximum allowable eccentricity e*. The line of apsides is rotated 
in such a manner that the solar pressure will cause the eccentricity to decrease to zero 
before increasing again. Figure 19 gives curves of e/fcp, w, and X when station 
keeping is not used. As seen from figure 19 the eccentricity is increasing when the 
apsidal line leads the Earth-Sun line (w - X > 0). The eccentricity is decreasing when 
the apsidal line lags the Earth- Sun line (w - X < 0). For method 3, the apsidal line is 
rotated so that a - x changes from the lead angle u>j - X^ to the lag angle u >2 - Xj. 

If the apsidal line is rotated through an angle Aw = 2(wj - Xj), at times t c and 3t c , 
then e/ep, on, and X curves as given in figure 20 are obtained. Appendix C shows that 
tangential thrusting requires only one-half as much AV to rotate the line of apsides as 
radial thrusting. For that reason, the requirements for this method are given for tan- 
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gential thrusting only. The AV per year for method 3 is 


Notice that 


AV 


s3 


pTT Vfl/l - /3^ 


3Skff 


2X \2 sinH/\ S in -1 /3 
2 > 


AV, 


s3 


= f^l - ? 


AV 


s2 


The acceleration level is 



( 35 ) 


(36) 


Figure 21 presents a plot of AV g g/k as a function of p with 0 as a family parameter. 
Figure 22 presents the corresponding plot for (M/k)a g g. 

Method 4 . - In method 4, the orientation of the apsidal line with respect to the 
Earth-Sun line is controlled such that the eccentricity is maintained at or slightly below 
the maximum allowable eccentricity e*. Figure 23 presents plots of e/ep, w, and X 
when method 4 is used. Initially the eccentricity is e*, and the apsidal line lags the 
Earth-Sun line by a small amount (w - X slightly less than 0). Due to the lag, eccen- 
tricity will decrease slightly. When the apsidal line coincides with the Earth-Sun line 
(u> = x), the eccentricity reaches a minimum. When the apsidal line leads the Earth-Sun 
line (w - X > 0), eccentricity increases and eventually becomes equal to e* again. At 
this time, the apsidal line is rotated so that u> - X has the same lag value it had at time 
equal to zero. The process is then repeated. Each time the eccentricity becomes equal 
to e*, the apsidal line is rotated through the same angle Acu. 

For very frequent corrections, the angle Ao> is very small. As a result, w « X 
and e « e*. Frequent corrections are desirable because the AV per year decreases 
as the frequency of correction increases. 

The requirements for this method are given for tangential thrusting only. The AV 
per year for frequent corrections is 
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(1 - 2/3) 


(37) 



If corrections are made every N days, then the acceleration level is 


(38) 


Figure 24 presents a plot of AV g ^/k as a function of p with /3 as a family parameter. 

Figures 25(a) and (b) present the corresponding plot for (M/k)a g ^ with N = 7 and N = 30. 

Comparison of methods . - To explicitly compare the different methods of controlling 

eccentricity, figure 26 plots AV ./k (for i = 2, 3, and 4) as a function of /3 with a duty 

cycle of 0. 01. When /3 is 0, AV_.A is 106 kilograms per meter per second (m/sec)/ 

2 

(m /kg) for methods 2 to 4. Equation (29) shows that AV g j/k is independent of /3 and 
has a value of 142 kilograms per meter per second. 

Several factors must be considered in choosing one of these methods of controlling 
eccentricity. If a cold gas system is used for station keeping, it is desirable to choose 
a method having a small AV per year. If /3 is large, method 4 requires significantly 
less AV per year (see fig. . 26). However, this method requires more complex or more 
frequent station-keeping maneuvers than methods 2 or 3. If ^ is close to 0, then the 
AV per year is nearly the same for methods 2 to 4. In this case, method 2 might be 
chosen because it is the simplest. If a low-thrust, high-specific-impulse system is 
used for station keeping, then it is not so critical to minimize the AV per year. In this 
case, either method 1 or 2 may be best. 



Interactions Between Station-Keeping Methods 

To this point, the station-keeping requirements for correcting effects due to solar 
pressure and triaxiality have been discussed separately. Since longitude excursions 
are due to both solar pressure and triaxiality, two allowable longitude excursions AL g 
and AL t can be specified in such a way that 

AL g + AL t = AL (39) 
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where AL is the maximum allowable longitude excursion due to both solar pressure 

and triaxiality. How AL g and AL^ are chosen depends on the methods of eccentricity 

control and triaxiality control to be used. Method 1 for controlling eccentricity is a 

special case because the eccentricity is always zero. AL is not an applicable param- 

s 

eter. Therefore, AL t may be set equal to any number less than or equal to AL. When 
using method 2, 3, or 4, the veolicity increment AV ., for i = 2, 3, or 4, is dependent 

S X 

on AL through the parameter 0. For each of these three methods, AV . becomes 
s si 

smaller as AL g becomes larger. On the other hand, AV^ is independent of AL t> A 
general rule can thus be stated that AV requirements are minimized by choosing AL 

s 

only slightly smaller than AL. 

To demonstrate the interaction of station keeping for triaxiality and station keeping 
for solar pressure, assume that method 2w?i chosen for controlling eccentricity. 

The spacecraft Is assumed to have an area- to- mass ratio of 0.15 square meter per 
kilogram and a reflectivity of 0. 3, so that k = 0. 195 square meter per kilogram. 

„ Choose y Q to be 45°. For AL = 0.30°, choose AL g =0.26° and AL t =0.04°. With 
these assumed parameters, triaxiality corrections would be made every 20 days, and 
eccentricity corrections would be made every 65 days. Figure 27 presents a plot of 
Ay* as a function of time when radial thrust is used for eccentricity control. The solid 

I L 

; line represents A y t . In any one orbit, Ay (in deg) is given approximately by 
£ Ay = Ay t + Ay g = Ay t + 2(57. 3e)sin 0 £ t (40) 

| where e is the instantaneous eccentricity of the orbit. Thus, Ay oscillates between 
£ Ayj. - 2(57. 3e) and Ay^. + 2(57. 3e). The dashed lines in figure 27 represent Ay^+2(57.3e) 
. and Ay t - 2(57. 3e). For example, on the 50 m orbit (t = 50 in fig. 27), Ay oscillates 
* between -0. 16° and +0. 24°. 

Tangential thrust for eccentricity control is more desirable than radial thrust be- 
cause only half as much AV is required (see eqs. (D2) and (D5)). However, unlike the 
radial-thrust case, tangential thrust for circularizing the orbit causes a change of 
3ire/4 radians in Ay t (see appendix D). Figure 28 presents a plot of Ay t as a function of 
time when tangential thrust is used. At t = 62, the first orbit circularization is made, 
causing an increase of 0.31° in Ay t . The triaxiality effect is used to advantage here 
because it causes Ay t to decrease steadily until Ay^ = -0. 04° on the 82nd orbit, at 
& which time the triaxiality correction is made, it should be pointed out that figures 27 
? and 28 are idealized curves used only to give a feel for how triaxiality, solar pressure, 

•y and east-west station keeping affect Ay. In an actual situation, other factors must be 
considered which might alter these curves. In particular, one must consider the effect 
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i 
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TABLE n. - STATION-KEEPING METHODS AND REQUIREMENTS 


Perturbation 


Requirements 


Sun and Moon Change i every 
N days 


Solar pressure 1 - Continuous thrust 


aj— -<2.3xl0-°)| 


/3SkA/4 


/ N/M \ 

r4) 


a sl ' Sk 


2 - Circularize whenever av^ = ( 


P* \/ P 


|2 sin Hly in~ L fy 


E\(A 


V M A 8X /^inH 


3 - Change whenever AV_ 


: /3Sk A/ pir \/p Vl - 

' ^ 'bsinnA Sin - 1 


SkV Jg E\/<3 f 1 ■ 


A sinH 


4 - Change w slightly 
whenever e = e* 


2 sin H, 

V ay 


. /SkV3N\ /l - a) 

\mA 4 /y 8ln H 


Triviality Change a whenever 
|Ay t | = AT r 


AV t /— i = 1.75 1 sin 2v I 
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of errors in north-south corrections. Misalinement of the north-south thrust vector may 
cause an in-plane error. From reference 2, the AV per year to correct these in-plane 
errors is of the order of three meters per second. 

The station keeping methods and requirements discussed in this section are sum- 
marized in table II. The equations given in table II are available in curve form in the 
report. Their use is best demonstrated by considering the sample problem given in 
appendix F. 


CONCLUDING REMARKS 

Station-keeping requirements have been determined for a class of high-power syn- 
chronous equatorial communication satellites characterized by large Sun-tracking solar 
arrays. The requirements for north-south control and for east-west control due to the 
Earth's triaxiality are the same as for previous communication satellites. However, be- 
cause of the larger solar arrays, the effect of solar radiation pressure must also be con- 
sidered for this new class of satellites. Solar radiation pressure produces an accelera- 
tion proportional to the area-to-mass ratio of the satellite, and the resultant accelera- 
tions change the eccentricity of the satellite orbit. The eccentricity causes an apparent 
daily east-west oscillation in the position of the satellite. For high area-to-mass ratio 
satellites, the east-west drift would be nearly i° longitude which would require station 
keeping for some missions. 

Equations and curves are given based on the assumption that all station-keeping 
corrections are carried out over a specified number of consecutive orbits with two 
rusting periods per orbit. For north-south corrections, the thrust is directed toward 
the north in one thrusting period and directed toward the south in the other. Alterna- 
tively, north-south corrections could be made with only one thrusting period per orbit, 
the thrust always being in the same direction. The curves giving AV for north-south 
corrections can be modified to handle the case of one thrusting period per orbit. The 
AV curves for east -west corrections, however, s ~e usable only for the case of two 
thrusting periods er orbit. 

This report covers the station keeping problem with emphasis on nonimpulsive 
low-thrust methods of station keeping. Parametric equations and curves giving AV and 
thruster acceleration requirements for the various methods of station keeping as a 
function of duty cycle and frequency of correction are presented. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 14, 1970, 

164-21. 
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APPENDIX A 


RADIO-FREQUENCY RADIATION PRESSURE 

The momentum of a quantum of energy is given by 

H = — = — 

C C 


(Al) 


where E is the energy of the quantum, C is the velocity of light, h is Planck’s con- 
stant and v is the frequency of radiation. If a plane source is emitting electromagnetic 
radiation, the force F associated with the radiation is given by 


F = 


A (H ) = JLM = w 

dt dt \C/ C 


(A2) 


where W is the total radiated power. The acceleration of a satellite due to radio- 
frequency radiation pressure is given by 


a 


F _ W/m 



(A3) 


Based on a solar array packing density of 100 watts per square meter, a maximum area- 
to-mass ratio of 0. 15 square meter per kilogram, and assuming that all of the collected 
power is radiated as radio-frequency power, then the maximum power- to- mass ratio 
W/m = 15 watts per kilogram. Hence, 

/ \ ^ Z m _ W = .A5 watts/kg 5x10 -8 m/sec 2 (M) 

' W C 3X10 8 m/sec 

A focused beam of radio-frequency energy emanating from an equatorial synchronous 
satellite may be deliberately pointed off the local vertical by as much as 8. 7°. The re- 
sultant acceleration vector a r j of such a beam may be resolved into in-plane and out-of- 
plane components. The out-of-plane component, being nearly constant, will have a neg- 
ligible effect on the orbit. The in-plane component of a rf may be further resolved into 
radial and tangential components. 

The radial component of a^, has the same kind of effect as the Earth's oblateness. 

If an adjustment in the orbit radius is made, no station keeping is necessary. The ad- 
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justment in radius for oblateness effects was 520 meters. For a radial component of 

__ _o 

a r j equal to 5x10” meter per second squared, the adjustment in radius for radio- 
frequency effects is only 3 meters. 

The effect of the tangential component of a rf is similar to that of triaxiality in that 

a steady longitudinal drift of the satellite is produced. Let a. denote the tangential 

n 

component of a_ f . The calculation of AV due to a. is the same as the calculation 
rt l rf 

of AV due to triaxiality. It can be shown that the AV is proportioral to the magnitude 

of the perturbing acceleration. The maximum value of a. would occur if the satellite 

n 

antenna is directed toward the east or west horizon at the equator. For this case, the 

_Q 

maximum value of a. is 7. 5x10 meter per second squared. The maximum tangen- 

rf 

tial acceleration due to triaxiality is seven times greater. The AV due to a f is no 

1 rf 

more than one-seventh as much as the AV due to worst-case triaxiality. The total AV 
per year due to radio-frequency effects does not exceed 0. 3 meter per second. 


23 



APPENDIX B 


ANALYSIS OF e(t) AND aMWHEN 3 Q = 0 

In this appendix equations are derived for e(t) and w(t) when the initial orbit is 
circular (e Q = 0). As mentioned previously, solar radiation pressure is the only pertur- 
bation having an appreciable effect on e and w. From reference 6, the time rate of 
change of e and u> is given by 


de _ -_3|^a(l -_e1 gk(Q • U) 

(Bl) 

dt 2y/J[ 


= sYa(l - e 2 ) a(? . fj, 

(B2) 

dt 2y/]Ie 


For small e, the equations can be simplified to yield 


de = -3Sk ( q . u) 

(B3) 

dt 2V 


dw = 3Sk (P . U) 

(B4) 

dt 2Ve 



To simplify the analysis, assume the orbit plane, equatorial plane, and ecliptic plane 
are one and the same. Then, as shown in figure 29, the unit vectors U, P, Q lie in 
this plane. Further, assume that the Sun is initially along the X-axis so that X = Xt. 
Expanding the dot products in equations (B3) and (B4) gives 

^ ^ sin(w - Xt) (B5) 

dt 2V 

4" = 3Sk cos(u , . ^ t) (B6) 

dt 2Ve 

Figure 30 shows the procet jy which solar pressure causes an initially circular orbit 
to become eccentric. When the satellite is in the vicinity of the positive Y-axis, solar 
pressure accelerates its motion and causes it to seek a higher altitude. When the satel- 
lite is in the vicinity of the negative Y-axis, solar pressure decelerates the satellite's 



% 


\ 


motion and causes it to seek a lower altitude. The result is a shif.ing of the orbit, 
creating a perigee on the positive Y-axis and an apogee on the negative Y-axis. There- 
fore, if e Q is 0, then u; o must be jt/ 2 radians. The solution to equations (B5) 
and (B6) with the initial conditions e Q = 0, o> 0 = vi/2 is found by assuming a linear vari- 
ation of a> with time. The solution is 


e = 


At 

2 


sin ^ t 

(B7) 

VA 2 


, 0 < t == 1 year 

(B8) 


/ 


\ 
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APPENDIX C 


ANALYSIS OP e(t) AND w(t)WHEN e Q * 0 

In appendix B, analytic solutions for e(t) and cu(t) were derived for the case e = 0. 

1 0 
Consider the initial conditions e Q = g ana - 0. Using the same assumptions and 

definitions as in appendix B, analytic solutions for e(t) and o>(t) can be found. The 

equations for e and w are 


e = (Q . u) ( C1 ) 

2V 


w = (P • U) (C2) 

2Ve 


Assuming that the vectors U, P, and Q lie in the same plane allows equations (Cl) 
and (C2) to be written as 


e = sia(ti) - At) 
2V 


(C3) 


u> = cos ( w _ * t ) (C4) 

2Ve 

The solution to equations (C3) and (C4) with the initial conditions e = \ e and w = 0 

o ^ p o 

is found by assuming a linear variation of a> with time. The solution is 


r _ 3Sk _ 1 p 
” 2VX ~ 2 P 


(C5) 


(0= Xt (C6) 

This solution for e(t) and w(t) corresponds to the Sun-oriented orbit. When e has a 

l o 

value other than 0 or ^ e p , equations (C3) and (C4) are not amenable to closed-form 

solutions. Computer solutions were obtained for e(t) and w(t) by numerically integrat- 
ing equations (Cl) and (C2). The assumption of the planar problem was not used when 
obtaining computer solutions, hi all cases, it was assumed that X Q = 0 (starting at 
autumnal equinox). For any values of e Q and o> 0 , e(t) and u>(t) were found to be peri- 
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odic functions with period of 1 year. In all cases, the apsidai line made either 0 or 1 net 
revolution per year. Figure 31 presents a plot of e/k as a function of time for the ini- 


tial conditions e = 0. 1 e and u> = 0. Notice that the minimum eccentricity is greater 
o P o 

than zero and the maximum eccentricity is less than e . For a given e Q , define two 
functions of Let e max(o> o ) be the maximum value of eccentricity obtained when 

(e Q , u> Q ) are the initial conditions, and let e min,a> 0 ) be the minimum value of eccen- 
tricity when (e Q , o> o ) are the initial conditions. For e Q = 0. 1 ep, figure 31 shows that 
e min(0) = 0. 002 k and e max(0) = 0. 018 k. Figures 32 to 34 present e max(w o )/k and 


e min(ci> o )/k for three different values of e . hi all cases, e max is smallest when 
u> 0 = 0. Of particular interest is the case e Q = ^ e (fig. 34). e max(0) and e min(0) 
are nearly equal, implying that the eccentricity is nearly constant. The reason for this 


is that the initial conditions 


e o = 2 e p’ 


co o = 0 


correspond to the Sun- oriented orbit. 





APPENDIX D 


STATION-KEEPING TECHNIQUES FOR CHANGING ORBITAL 
ECCENTRICITY AND ROTATING THE LINE OF APSIDES 


Circularizing a Slightly Eccentric Orbit 

In this appendix, it is assumed that all station-keeping corrections are completed 

in a 24-hour period (M = 1). A slightly eccentric orbit can be circularized by radial 

thrusting (thrusting in a direction collinear with the radius vector) or by tangential 

thrusting (thrusting in a direction collinear with the orbit velocity vector). Figure 35 

presents a sketch of the radial thrusting maneuver when impulsive thrusting is used. 

Position Aj corresponds to a true anomaly of 90°. This correction is made by first 

thrusting inward (toward the Earth) in the vicinity of A, and then thrusting outward in 

1 1 

the vicinity of Ag. The transfer orbit has an eccentricity of ^ e. The AV C for impul- 
sive thrusting is 


AV = eV (Dl) 

c i 

When nonimpulsive thrusting is used, 

AV = eV — ^ — (D2) 

2sinH 

2 

Let y be the value of y at the perigee of the orbit just prior to the circularization 
maneuver. Then define Ay by the equation 

r 

Ay p = y - y p (D3) 

where y p is the mean longitude in the sense that the satellite longitude (before the 
station- keeping correction) oscillates about y p with an amplitude of 2e radians and a 
period of 24 hours. Figure 36 presents a plot of Ay p as a function of time when the 
duty cycle of the station-keeping correction is 0.01. Notice that Ay p is zero after the 
completion of the station-keeping maneuver. 

Figure 37 presents a sketch of the tangential thrusting maneuver when impulsive 
thrusting is used. Position Aj is the orbit perigee for tangential thrusting. The cor- 
rection is made by first thrusting westward in the vicinity of A ^ and then thrusting 
eastward in the vicinity of Ag. The transfer orbit has eccentricity ^ e. 


/ 


1 
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for impulsive thrusting is 


The AV 

c 


AV 


c 



(D4) 


When nonimpulsive thrusting is used, 

AV = - eV — 21 — (D5) 

c 2 0 . p7T 

2 sin *— 

2 

which is one-half as much AV c as required with the radial correction scheme. 

Figure 38 presents a plot of Ay^ as a function of time when tangential thrusting is 
used. Notice that a residual longitude error of 3rrt/4 radians i« left after the comple- 
tion of this station-keeping maneuver. It can be shown that an additional AV of approxi- 
mately eV/4d is required to remove the residual longitude error if d is the time in 
days to make the correction. 


Rotating Line of Apsides 

Now consider radial and tangential thrust maneuvers which will rotate the line of 
apsides of an orbit having a small eccentricity e. Let \I/ be the angle through which 
the apsidal line is to be rotated. Figure 39 is a sketch of a radial correction maneuver 
when impulsive thrusting is used. The correction is made by thrusting outward in the 
vicinity of and thrusting inward in the vicinity of Ag. With radial thrusting, Aj 
corresponds to a true anomaly of ^ ^ in the final orbit. When impulsive thrusting is 
used, AV C is given by 


AV , = 2eV sin ^ 
L 2 


For nonimpulsive thrusting, AV C is given by 


AV = 2e\ sin^f-J* 

C 2 \2sin Hj 

2/ 


(D6) 


(D7) 
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Figure 40 presents a plot of Ay^ as a function of time when rotating the apsidal line 
20° with radial thrusting. Notice that Ay^ is zero after the completion of the station- 
keeping maneuver. 

When rotating the line of apsides with tangential thrusting, corresponds to a true 
anomaly of -g + 90° in the final orbit. The correction is made by first thrusting west- 
ward in the vicinity of Aj and then thrusting eastward in the vicinity of Ag. Figure 41 
is a sketch of this maneuver when impulsive thrusting is used. 

When impulsive thrusting is used, AV C is given by 

AV = eV sin ^ (D8) 

c 2 


For nonimpulsive thrusting, AV C is given by 


AV = eV sin ^1 — E!_ 

2 \2sinH 
2 


(D9) 


Only half as much AV„ is required with tangential thrusting as with radial thrusting. 

V * 

Figure 42 is a plot of Ay p as a function of time when rotating the apsidal line 20° with 
tangential thrusting. A residual longitude error of (3ne/2) sin ^ radians is left after 
the completion of the maneuver. 
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APPENDIX E 


DERIVATION OF EQUATIONS FOR AV sj AND a sj FOR THE 
FOUR METHODS OF CONTROLLING ECCENTRICITY 

METHOD 1 


Method 1 is to continuously cancel the effect of solar pressure by thrusting in a di 
rection toward the Sun. The acceleration level is given simply by 


a gl = Sk (Ei) 

The AV per year is given by 



It will be found convenient to express AV g ^ as 



METHOD 2 


Method 2 is to circularize the orbit each time the eccentricity becomes equal to e*. 
For this method, AV (AV per correction) can be found by using equation (D5). 



Recalling that e* = j3e p and from equation (C5) 

e .» 

P VX 


(E4) 
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allows equation (E4) to be written as 




\ 


AV - 3S ^ ( 

C \2sinHj 

2/ 


From appendix B equation (B7), the unkept eccentricity is 


e = e„ sin - t 

P 2 


The time t c between corrections can then be found from the equation 


e_j3 = e* = e„ sin - t 
P P 2 ' 


Solving for t yields 


t _ 2 sin' 1 ^ 

c " x 


The number of corrections per year K is given by 


_ 1 year _ (2 
t n " \> 


2tA 1 12 sin" 1 ^ jr 


sin’ 1 ^ 


AV g 2 is then given by (K)(AV C ). 



A 


2 sin P^y\sin~*/3/ 


The acceleration level can be found from the equation 


av c - saF™ 


% 
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Solving for a g p yields 



(E12) 


METHOD 3 


Method 3 is to rotate the line of apsides in such a manner that the solar pressure 
will cause the eccentricity to decrease to zero before increasing again. Let 2t c be the 
time between corrections. Using the same derivation as in method 2, 


. 2 sin _1 /3 

c " X 

AV C can be found by using equation (D9) to be 



The values of w and X at time t are 

c 

w = — t + — = sin"*/3 + - 
2 c 2 2 

X = Xt = 2 sin _1 /3 

u 


The angle Aw is then given by 


Aw = 2 (w - X) = jt - 2 sin" 


(E13) 


(E14) 


(E15) 

(E16) 


(E17) 
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Solving for a 


gives 



(E23) 


METHOD 4 

Method 4 maintains the eccentricity nearly equal to e* by frequently rotating the 
apsidal line in such a manner that co « X. To derive formulas for and a g ^, 

return to the planar problem where it was assumed that the Earth, the Sun, and the sat- 
ellite are in the same plane. Further assume that the Sun and the orbit perigee initially 
are along the X-axis and that /3 < -g. Thus, X = Xt and u> 0 = 0. In equations (B 5) 
and (B6), e and cl> are given by 


sin(co - Xt) 

(£24) 

2V 


cos(co - Xt) 

(E25) 


2Ve 


With method 4, the apsidal line is rotated a small amount each time the eccentricity 
becomes eqv to e*. This small rotation would cause the eccentricity t j decrease only 
slightly b* > increasing again. Eccentricity as a function of time for this method is 
shown in figure 43. In this curve, 2f is the time between corrections, and ae* is the 
minimum eccentricity. The time t will be relatively small, and a will be only 
slightly less than 1. Since e is kept nearly constant, then de/dt must be approxi- 
mately zero. From equation (E24), it follows that must be kept nearly equal to Xt. 
Thus the apsidal line must be controlled so that perigee remains directed toward the 
Sun. Assume now that at t = 0, e = ae*. From equation (E25), for small t, 


• 3Sk 

co = 

2Ve „ 


(E26) 


By assumption e Q < -j e^. Using the equation ‘j e p = 3Sk/2VX, we obtain w > X. The 
apsidal line, if uncorrected, will rotate faster than the Sun. The control of the apsidal 
line must then be as given in figure 43. 
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The times of correction are t , 3t 5t , . . . . With the assumption that t is 

c c c c 

small, AV requirements can be determined by linearizing and solving equations (E24) 
and (E25) in the time interval 0 i t s t . The initial conditions are =0, X Q = 0, and 
e Q = are * . The linearized equations are 

^ = 3Sk (w _ ^ t) (E27) 

dt 2V 


dcu _ 3Sk 
dt 2Ve Q 


(F28) 


Equation (E28) can be further simplified by substituting for e Q : 

e _ 3Sk ag ( E 29) 

0 VX 

U.ing the assumption that or is approximately equal to 1, equation (E28) becomes 

^ = A. (E30) 

dt 2/3 


Solving for u>, 



Substituting this solution for o> into equation (E27) gives 

de _ 3SkX /l _ _ A t 
dt 2V \2/3 ) 


which can be integrated to give 


e 


e o + 


askx l 


,2/3 



(E31) 


(E32) 


4V 


(E33) 



I 


From equation (E33), t is calculated to be 


= 2 '£ / 2(1 - a) 
c ' X V 1 - 2 p 


(E34) 


The rotation angle Aw of the apsidal line is 


Aw = 2(wt - Xt) = 2xt (LlJ£ )= 2 V2(l - o)(l 
c c c \ 2/3 / 


2 / 3 ) 


(E35) 


Assuming Aw is small and a is approximately 1, AV c can be found by using equa 
tion (D9) 


AV 


c 



AV 


c 



Rearranging terms, 


AV c 


3Skt^ 

2 



The number of corrections per year K is given by 

k = 1 y ear = JL 

2t Xt 
c c 


(E36) 


(E37) 


(E38) 


(E39) 


AV g4 is then given by (K)(AV C ). 
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(E40) 


(E41) 


(E42) 


(E43) 
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APPENDIX F 


SAMPLE PROBLEM 
INTRODUCTION 

To demonstrate the use of the curves in this report an example of a typical high- 


power communication satellite mission is presented. ne satellite parameters are 

Mass, kg 1000 

2 

Area-to-mass ratio, m /kg , 0.154 

Average reflectivity , 0.3 

The mission requirements are 

Allowable longitudinal error, deg 0.2 

Longitude position, deg west longitude 95 

Mission life, yr (does not include 1 yr reserve) 1.5 

Thruster systems available are 

(1) Thrust, N 9. 8 

Specific impulse, I , sec 100 

(2) Thrust, N 4. 9x10' 3 

Specific impulse, I , sec 2000 


-3 

The high-thrust system yields an acceleration of 10 g on the satellite; the low- 
thrust system yields an acceleration of 5x10" g on the satellite. 

NORTH-SOUTH STATION KEEPING 
High -Thrust System 

Figures 11 and 12 are used to determine AV^g. The firnt step is to determine 

from figure 12 the duty cycle p for given values of M, N, and a^g. Assume that 

corrections are made once every 60 days (N = 60) and that the correction is carried out 

— 3 

in 1 day (M = 1). The ordinate Ma^g of figure 12 is 10 g, and the corresponding 
abscissa p is approximately 0. 01. So the correction is accomplished by thrusting in 
the northerly direction for approximately 0, 12 hour, and a half orbit later thrusting in 
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the southerly direction for approximately 0. 12 hour. Knowing p from figure 12, the 
AV per year, AVj^g, is then determined from figure 11 to be 46 meters per second. 


Low-Thrust System 

For the low-thrust case, assume first that corrections are made daily (M = 1, 

_ ^ 

N = 1). The ordinate, Ma^g, of figure 12 is 5x10 g, and the corresponding abscissa p 
is 0. 3. From figure 11, AV N g is 48 meters per second. If it is now assumed that cor- 
rections are made weekly (N = 7), the smallest possible value of the ordinate Ma^g in 
figure 12 is approximately 1. 5xl0~ 6 g, corresponding to a unity duty cycle p. Since 

_7 

a^g is 5x10 g, M must be 3 for this case. Thus the correction scheme consists of 
continuous thrusting for 3 days and no thrusting for the next 4 days. Since the duty cycle 
p is unity, the AV per year from figure 11 is 72 meters per second. 

Figures 11 and 12 can be modified to handle the case of one thrusting period per 
orbit instead of two. The modifications are: 

(1) The thrusting time per orbit divided by the orbit period is 2p (instead of p). 

(2) The number of days between the beginnings of successive inclination corrections 
is N/2 (instead of N). 


EAST-WEST STATION KEEPING DUE TO TRIAXIALITY 

The AV required for triaxiality is a function only of y Q . For a desired satellite 
position of 95° west longitude, y Q = 105° - 95° = 10°. The AV per year is given in 
figure 15 to be 0. 5 meters per second. 


EAST-WEST STATION KEEPING DUE TO SOLAR PRESSURE 

Method 1 

For continuous thrusting against the Sun, the AV per year, AV g j, is given by 
equation (29) to be 

AV »>‘ir(s) (F,) 

For this satellite, the AV per year is 28. 4 meters per second. 
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Method 2 


High-thrust system. - In determining AV requirements lor the remaining station- 
keeping methods, AL g and ALj. must first be chosen. Let AL g = 0.15° and 
ALj. = 0. 05°. The parameters k and fi can be calculated from equations (4) and (32), 
respectively. 


<3 = 0. 4 


> A =0.20 HL 

m kg 

(F2) 

AL s 

2 = 0. 30 

(F3) 


To determine AV g2 , refer to figures 17 and 18. The first step is to determine from 
figure 18 the duty cycle p for given values of M and k. Assuming M - 1, the ordi- 
nate (M/k) a g 2 of figure 18 is 5x10 g's per square meter per kilogram. The corre- 
sponding abscissa p, with the family parameter /3 = 0. 30, is seen to be much less than 
0. 01. Thus the orbit is circularized by two thrust pulses (12 hr apart) of duration much 
less than 0. 24 hour. Knowing p and fi, figure 17 shows AV^/k is 104 (m/sec)/(m^/kg), 
so that the AV per year, AV s2 , is 20. 8 meters per second. 

Low-thrust system . - From figure 18, with the family parameter /3 = 0. 30, the 
smallest possible value of the ordinate (M/k)a g 2 is approximately 2. 0x10"® g's per 
square meter per kilogram, corresponding to an abscissa p of 1. 0. Since (l/k)a g 2 is 
2. 5xl0” 6 g’s per square meter per kilogram, M is 8 in this particular case. Thus, for 
this low-thrust case, 8 days of continuous thrust is required to circularize the orbit. 
From equation (E8), corrections are made once every 35 days. For a unity duty cycle 
and /3 = 0. 30, figure 17 shows AV g2 A is 162 (m/sec)/(m^/kg), so that the AV per 
year, AV g2 , is 32. 4 meters per second. 

Since AV s2 is largest when p=1.0, it is desirable to choose a smaller p. If p 
is chosen to be 0. 3, then from figure 18, (M/k)a g2 is 3x10"® g's per square meter per 
kilogram, implying that M is 12. For this case, the orbit is circularized by thrusting 
7. 2 hours per orbit for 12 consecutive orbits. From figure 17, the AV per year, 

AV g2 , is 21. 2 meters per second, considerably less than the 32. 4 meters per second 
for a unity duty cycle. 
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Method 3 


High-thrust system . - To determine AV g;j , refer to figures 21 and 22. Assuming 
M = 1, the ordinate (M/k)a s g of figure 22 is 5x10” ^ g's per square meter per kilogram. 
The corresponding abscissa p, with the family parameter 0 = 0. 30, is seen to be less 
than 0.01. Thus the apsidal line is rotated by two thrust pulses (12 hr apart) of duration 
less than 0. 24 hours. Knowing p and 0, figure 21 shows AV g gA is 100 (m/sec)/ 
(m 2 /kg)> so that the AV per year, AV g g, is 20. 0 meters per second. 

Low-thrust system . - If p is chosen to be 0. 4, then from figure 22 (M/k)a g g is 
6 x 10“5 g’s per square meter per kilogram, implying that M is 24. For this case, the 
apsidal line is rotated by thrusting 9. 6 hours per orbit for 24 consecutive orbits. From 
equation (E13), corrections are made once every 70 days. From figure 21, AV s g/k" is 
105 (m/sec)/(m 2 /kg), so that the AV p°r year, AV^^is 21. 0 meters per second. 


Method 4 

High-thrust system . - Unlike methods 2 and 3, the parameter N in method 4 is not 
a function of 0. The only restriction on N in method 4 is that it be small enough to 
justify linearising the differential equations for e and w. For the range of parameters 
considered, the linearized equations can be justified for N £ 30. To determine AV g4 
for the high-thrust case, refer to figures 24 and 25. Assuming N = 30 and M = 1, the 

_ O 

ordinate (M/k)a g ^ of figure 25(b) is 5x10 g’s per square meter per kilogram. The 
corresponding abscissa u, w\th the family parameter 0 = 0. 30, is seen to be much less 
than 0.01. Thus the apsidal line is rotated by two thrust pulses (12 hr apart) of duration 
much less than 0. 24 hour. Knowing p and 0, figure 24 shows AV g4 A is 42 (m/sec)/ 
(m 2 Ag), so that the AV per year, AV g4 , is 8. 4 meters per second. 

Low-thrust system. - If N is now chosen to be 7, and a duty cycle p of 0.2 is 
desired, then from figure 25(a), (MA)a g4 is approximately 5xl0~° g’s per square 
meter per kilogram, implying that M is 2. For this case, the apsidal line is rotated 
by thrusting 4. 8 hours per orbit for 2 consecutive orbits. From figure 24, AV g4 A is 
43 (m/sec)/(m 2 Ag), so that the AV per year, AV g4 , is 8.6 meters per second. 

SUMMARY OF REQUIREMENTS 

Assuming the ratio of propellant mass to spacecraft mass is small, the propellant 
mass as function of AV is given by 
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% 


m p 


m AV 



(F4) 


where m^ is the propellant mass, m is the spacec art mass, and g is the accelera- 
tion of gravity (9. 8 m/sec 2 ). Equation (F4) can be used to calculate the propellant mass 
once the AV is known. Table m summarizes the station-keeping requirements. 


TABLE ffl. - SAMPLE PROBLEM STATION-KEEPING REQUIREMENTS FOR 5 + 1 YEAR MISSION 



High thrust 

Low thrust | 

Propellant 

mass, 

kg 

Velocity 

increment, 

AV, 

m/sec 

Duty 

cycle, 

P 

Time 

between 

corrections, 

days 

Propellant 

mass, 

kg 

Velocity 

increment, 

AV, 

m/sec 

Duty 

cycle, 

P 

Time 

between 

corrections, 

days 

North- south 

282 

276 

0.01 

GO 

14.7 

288 

0.30 

1 

Solar pressure 
method 4 

51 

50 

<0.01 

30 

2.7 

52 

0.20 

7 

Triaxiality 

3 

3 



°' 2 ] 

3 

— 
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Figure 4, - Variation in satellite longitude as function of time. Perturbation, triaxiality; desired sat 
ellite longitude measured from Earth's minor axis, 45°. 
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Figure 5. - Inertial coordinate system. 
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Figure 6 . - Normalized eccentricity as a function of time. Initial conditions: eccen- 
tricity, a 



Figure 7. - Longitude of perigee and longitude of Sun a? functions of time. Initial 
conditions: eccentricity. 0; longitude of Sun, a 
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Figure 8. - Normalized eccentricity as a function of time. Initial conditions.* eccen- 
tricity, 1/2 e p ; longitude of perigee, CL 



Figure 9. - Longitude of perigee and longitude of Sun as functions of time. Initial 
conditions: eccentricity, 1/2 e p ; longitude of perigee, C ; longitude of Sun, a 
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Figure 24. - Normalized velocity increment per year for station keeping with method 4 
as function of duty cycle with eccentricity ratio as a parameter. 
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Figure 26. - Normalized velocity increments per year for methods 2 , 3, and 4 as func 
tions of eccentricity ratio. Duty cycle, a 01 



Figure 27. - Variation in satellite longitude as function of time when 
station -keeping method 2 with radial thrusting is used. 


60 













Initial perigee longitude, u) Q , deg 

Figure 34. - Maximum and minimum normalized eccentricity as function of initial perigee 
longitude. Initial conditions: eccentricity, CL5e p ; longitude of Sun, a 
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Figure 36. - Variation in satellite longitude as function of time when circularizing orbit with radial 
thrust Eccemricity of original orbit is e. 
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Figure 37. - Circularizing the orbit using two tangential impulses. 
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